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1. INTRODUCTION 

Wind energy in addition to photovoltaic energy source is the fastest growing energy with an annual 
growth rate of 30% and with a predictable penetration of 12% of the global electricity demand [1]. The use of 
power electronic converters allows for a variable speed operation of the wind turbine where the wind energy 
conversion system extracts maximum power from the turbine during peak operation as reported in [2]. One 
of the energy generation systems commercially available in the wind energy market is the doubly fed 
induction generator (DFIG) and its numerous advantages are enumerated in [2], [3]. To produce a maximum 
energy, an efficient DFIG with higher power rating is a precondition for wind energy conversion system 
(WECS) though active power loss during the operation is the main drawback for efficiency optimization and 
life expectancy’s challenge. Variable speed wind turbines are usually more effective due to their improved 
efficiency in capturing more wind power and their innate ability to achieve higher quality of power at 
optimum wind speed [4]. DFIG implementation is increasing in leaps and bounds as a consequence of a 
reduced mechanical stress and noise in addition to the flexible control of active and reactive power which is 
based on the back-to-back power electronics converter sandwiched between the induction machine and the 
power grid [5]. Active power generation at the stator and rotor terminals is effectively controlled using the 
rotor side converter (RSC). The regulation of the stator side active and reactive power is independently 
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achieved with the RSC. However, the fluctuations in power output from the DFIG-based wind generation in 
response to the variations in wind speed adversely affect the needed power quality [6]. Therefore, with an 
increased penetration of the DFIG in power system operation, inertia of the power system is reduced. An 
increased capacity in grid integration of this variable power poses an impactful challenge on power system 
stability [7]. In line with the system stability is the minimization of the overall system losses under different 
wind speed and reactive power conditions. Therefore, to improve on the system overall efficiency, the 
reactive power flow is regulated with a minimized system electrical loss on the generator and on the power 
electronic converters. 

In this paper, a modified equivalent circuit of the DFIG with variable core resistance values for loss 
minimization was adopted. The copper and iron losses of the DFIG were modeled as a function of the rotor 
dq-axes currents and stator flux. At optimum condition, the stator d-axis flux that minimizes the total DFIG 
loss was derived. The simulation results obtained showed that the efficiency values of 83.45% and 41.21% 
were obtained when the equivalent paralleled core resistance values of the DFIG were kept at 0.75Q and 
0.25Q. The power loss at a slight transient disturbance peaked at 202.277 kW and reduced to 17.2058 kW 
under steady state for 0.75Q net resistance. Similarly, at a net core resistance value of 0.25Q, the power loss 
during a transient disturbance peaked at 757.325 kW and reduced to 23.7597 kW. The variations in the power 
loss at different net-core resistance accounted for the difference in the efficiency values as presented in the 
preceding simulation waveforms. 

This paper is organized as follows: Section 1 is the introduction; Section 2 is the reviewed literature 
on DFIG loss control; Section 3 illustrates the wind turbine model of doubly fed induction generator; 
Section 4 presents the DFIG mathematical modeling and power loss minimization scheme; and Section 5 
contains the simulation results. Section 6 is the conclusion. 


2. METHODOLOGY AND REVIEWED LITERATURE ON DFIG POWER LOSS CONTROL 

Research procedures and methods applied in this paper involved mathematical modeling of the doubly 
fed induction generator at variable core resistances to determine its effects on the overall power loss of the 
system. Stability test using the Sylvester criterion method and data acquisition derived from simulation 
parameters were parts of the methods applied in the empirical analysis of the system controllability. In [8]—[10], 
direct power control (DPC) and direct torque control of the DFIG was proposed as a conventional control 
scheme. However, large torque and power ripples are the two major draw backs associated with this method. 
The various methods involving back stepping direct power control strategy for DFIG have been presented in 
[11]-[14]. Although a detailed analysis was presented but emphasis on variable core resistance and minimum 
flux at which the DFIG power loss is minimized was not presented. 

In [15], [16], a sliding mode approach in minimizing copper loss was presented in a nonlinear 
control application while a coordinated predictive control of the DFIG using non-Gaussian wind power 
predictive distribution was reported in [17]. The reports concentrated on regulating the extracted power from 
the wind turbine while incorporating machine loss minimization technique for iron-loss. In [18], [19], 
minimization of the active and reactive power ripples using the direct power control of matrix converter fed 
DFIG was presented without reference to the DFIG core losses minimization. 

This paper therefore analyzed a complete loss minimization of the DFIG using the parameter 
variation that involves the core resistance and stator flux. A complete DFIG loss equation was modeled and 
derived from a modified equivalent circuit involving the variable core resistance while a minimal flux 
equation that minimized the DFIG losses was also derived. The Sylvester criterion using the Hessian matrix 
was applied to determine if the DFIG total power loss was minimized within the context of the variable core 
resistance values chosen. Therefore, with the leading principal determinants positive for Req = Rfe = 0.25Q 
and 0.75Q, it showed that the overall total loss of the DFIG is controllable and so is minimized. 


3. WIND TURBINE MODEL OF DOUBLY FED INDUCTION GENERATOR. 

The schematic diagram of the DFIG wind turbines operation is presented in Figure 1. The DFIG 
stator is directly connected to the grid through the grid side converter GSC and its rotor terminal is connected 
to the rotor side converter RSC through the slip rings. The wind turbine is modeled with reference to optimal 
power tracking to provide a maximum energy capture from the wind. The aerodynamic model of the wind 
turbine gives a coupling between the wind speed and the mechanical torque produced by the wind turbine. 
Aerodynamic is a science of physical laws that deals with the behavior of objects in airflow and forces that 
are produced by this airflow. The aerodynamic power equation extracted from the wind turbines is given by 
(1) as reported in [20]-[22]. 
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Where Pm (W) is the mechanical power of the turbine, p (Kg/m*) is the air density, A (mô) is the area 
covered by the rotor turbine, V (m/s) is the wind speed upstream of the rotor and C, is the performance 
coefficient or power coefficient. The power coefficient is a function of the pitch angle (B) of the rotor blades. 
The tip speed ratio is the ratio between the blade tip speed and wind speed upstream of the rotor. Altering the 
pitch angle implies slightly rotating the turbine blades along the horizontal or vertical axis. Cp which 
represents the wind turbine power coefficient is given by (2). 


12.5 


Cp, B) = 22 x = — 0.0048 — 0.05) eM (2) 
1 1 _ 0.035 
AU A+0.08ß B3 +1 (3) 


Where B represents the blade pitch angle and à represents the tip speed ratio which is given by (4). 


= ork 
A = oe (4) 


The dynamic model wind turbine is associated with the rotor speed w, and gear box ratio nz by (5). 


_ optimum 
W, = Ng X o (5) 


The exact dynamic model of the torque equation for the generator is given by (6). 


= 
B 

| 

l 


3 
= TE x V2 x GA) (6) 


Where p = air density (Kgm™); R = radius of the turbine (m); Vw= wind speed (MS"!); Ci (A) = Torque 
coefficient (pu). The wind turbine characteristics are presented in Figure 2. 

The plot of power coefficient against the tip speed ratio at varying blade pitch angle is represented 
in Figure 2. It is observed that the power coefficient increases as the pitch angle decreases. Therefore, at a 
zero-pitch angle, the maximum power coefficient is 0.5 which implies that the maximum power is tracked at 
a zero-blade pitch and tip speed ratio value of ten as shown in Figure 2. 


Stator P&Q Transformer 


Figure 1. Schematic diagram of a DFIG operated wind turbine 
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Modified voltage equations can be obtained by substituting (11)—(14) into (7)-(10). 
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Wind Turbine Power Coefficient Cp against tip-speed ratio at varying Pitch Angle 
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Figure 2. A plot of wind turbine power coefficient against tip-speed ratio at varying pitch angle 


DFIG MATHEMATICAL MODELING AND POWER LOSS MINIMIZATION SCHEME. 


Stator and rotor voltage as well as flux equations as reported in [23], [24] are given in (7)—-(14) 


Vas = Rsias 
Vas = Rgigs 
Var = Ryigr 
Var = Rriqgr 


Aas = Lisias 
Ags = Lysigs 


dar = Lyriar 


Agr => Lurigr + Linigm 


Vas = Rsias 
Vas = Rsiqs 
Var = Rriar 
Vor = Rrigr 
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mn WeAgs 
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Modified equivalent circuit model of the DFIG with multiple paralleled core resistance is obtained from 
(15)-(18) and presented in Figures 3 and 4. Branch currents are obtained from Figures 3 and 4 by applying 
KCL as expressed in (19) and (20). 


igs + iar = 


idm + 


[-weLmigqm] Z 


iam — Ai 
Req dm qm 


(19) 
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+ [weLmidam] 


igs + igr = İqm Reg = igm + Aigm (20) 
A simultaneous solution of (19) and (20) by substitution gives rise to (21) and (22). 

. [Gastiar) +A(igs R iqr)] 

iin TO 21) 

f [(igs tigr) +A(ids + iar)] 

i age (22) 
Where: A = 2°“. For vector controlled condition, Ag; = 0 and Aas = As substituting this into (11) and 


eq 


(12) gives rise to (23) and (24). 
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Figure 3. D-axis equivalent circuit of a DFIG with multiple paralleled core resistance 
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Figure 4. Q-axis equivalent circuit of a DFIG with multiple paralleled core resistance 


Therefore, the dq-axes magnetizing currents under vector control is obtained by substituting (23) 
and (24) into (21) and (22) as re-presented in (25) and (26). 


[As+(Lis—Lm)(iar + Aiqr)] 


iam 
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[(LLs-Lm)(iqr- Aidr)—Aàs] 


lqm = LLs(1 + A?) 
The dq-axes magnetizing voltages are presented in (27) and (28). 
= f __ AWweLmàÀ @eLmbis(Aigr—igr) _ f ; 
Vam = —®elmigm = iis tay t Geen ANTA (Aiar — iar) < 
= y _  Welmà: ®eLmLis (iar + Aigr) _ - : 
Vam = —®elinlam = nara) oP ian. ee (iar + Aigr) ae 


WeLm Ss 
Where: Ay = nets + azyendA2 = LysA1- 

The voltage magnitude across the magnetizing branch is presented in (29). 

VA = Vin + Vem = A3(1 + AD [ihs + i2-] + 2A,A2(1 + A?)Jàsiar + AZ(1 + A?)A2 (29) 


Total core loss in terms of magnetizing voltage and equivalent core resistance is given by (30). 


(Vm)? f . 
PRR = GR = RACI + APGE) + Ag (L + A*)Aslar + AFC + A?) (30) 


rotor copper losses with core losses which are the total resistive losses of the DFIG are presented in (31). 


Where: i? = iĝ, + iâr and Reg = Rre||Ree |RrellRre n+1 are derived from Figures 3 and 4. The stator and 


1 = : g g p ` ; 
Po tive = Rs lis + i2] + R, [iã, + i2] + Rreližfe + Brel (31) 


Substituting (23) and (24) into (31) gives rise to (32). 


Total _ Rs (As)? RsLin r, \2 _ 2RskLmy ; RsLin f. 2 +2 +2 +2 
P Resistive — L2 + LÊ Gar) x LŽ Às lar + L2 (iar) + R, ir + Rfe [ize + ide | (32) 
Ss S S sS 


Other types of power losses associated with the DFIG operated system are the RL-filter losses, 
Frictional losses and power electronics converter losses as presented in (33)—(35) respectively. 


= +2 +2 

PRL-Fiter = Reitter(id—riter + ij Filter) (33) 
= 2 

Prrictiion = Kpr Wmech (34) 
= i2 12 

Peonverter a Reonverter lld- converter + ie zconveriter) (35) 


The total losses of the DFIG scheme can be written in terms of the total core loss, resistive loss, RL-filter 
losses, Frictional losses and power electronics converter losses as given in (36). 

Total _ p Total Total = 

P Loss P dre +P Resistive + PRL-Filter + Pfrictiion + Pconverter = 


, Rs (As)? 
2A, AzAsigr(1 + A?) + A325 (1 + A”)] + ae 


1 x 
= — [42 (1 + Ai? + 
Req 
RsL¥, .2 2RsLm 4; 2 .2 .2 
L2 lr 7 L2 Às lar + R, ir + Ree lidte + iåfel + 
Ls Ls 


PA i2 2 «2 2 
Reiter (ij-Fitter + iĝ filter) + Kpr Wmech + Rëohverter laeeonvente + Ta converter) (36) 


The flux at which Total Power loss is minimized is determined by setting the derivative of (36) to zero. 


dP Toss 1 . 2ÀsRs  2RsLm . 
dh, =0= Reg [2A,Azig-(1 + A?) + 20,47 (1 + A?)] + T — iz, iar = 0 
= [RrReqlm—(1+A7)A1A2Lisliar 
Asmin = A2(1+A2)L?.,+ RsReq G2 


The flux at which minimum loss occurred is as presented in (37). This implies that the minimum flux 
increases with a proportionate rise in the rotor direct axis current and equivalent core resistance values. 
Therefore, the overall total loss is minimized when A smin is substituted for A, in (36) that gives rise to (38). 
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Total — pTotal Total =o 2 
P Loss-min 7 P core—min +P Resistive—min g3 PRL-Filter + Pfrictiion + Pconverter m Reg [A (1 + 
2) 32 4 2 232 2 Rs(a smin)? RsLin 22 2RsLm . 
A ) Ip + 2A,A2A sminigr (1 +A ) F AiÀSmin (1 +A )] F L2 F L2 lr 7 L2 Asmin lar + 
Ls Ls Ls 
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+2 
ig onyen) (38) 


The DFIG power loss minimization technique applied on (38) was actualized with the Sylvester’s 
criterion for positive definite state using the Hessian matrix determinant test. The objective function to be 
minimized is the overall total power loss already given in (38). The Hessian matrix is obtained by taking the 
second derivatives of (38) in terms of stator flux and rotor current as presented (39)-(42). 


d2 Total 


z PLoss—min a 4 2) 42 2Rs 
Hi a Se AEEA AE (39) 
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d2protal hin 1 212,Rs 
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sS 


H22 = 


The Hessian matrix equation is presented in (43) for the DFIG total loss parameter stability test. 


Hy, = 


2Rs 
[201 + A*)At] + = H2 = 0 
H = 


eae 

Re 

i 1 2L2.R (43) 
Hz, = 0 H22 = —[2(1 + A7)A3] + == + 2R; 


Req L2 


The total power loss is minimized if the leading principal determinants of (43) are positive and greater than 


zero. The principal determinants are obtained with the DFIG parameters in Table 2 substituted to give the 
following: 


1 2) 42 2Rs 
[Hii] = [20 + AAG] + 45] = 691.9801 

eq s 

1 2R; 
[2(1 + 4°A5i] + — 0 
Hı Aya] _ Req Ls 
H Hal 
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1 212.R 
[2(1 + A2)A3] + 3 = + 2R; 
eq s 


1 2R 212,R 
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eq s S 

= 2.9156e + 05 


The parameters substituted above showed that the DFIG power loss model was controllable with a 
positive definite value of 691.9801 and 2.9156 e ** for the leading principal determinants of the Hessian 
matrix. For dynamic simulation analysis of the doubly fed induction generator under a load disturbance, the 
state space model equations were derived with the aid of Figures 5 and 6 and presented in (44)—(57). 
Applying the Kirchhoff’s current law at nodes A and B on the simplified DFIG dq-axes equivalent circuit 
diagram presented in Figures 5 and 6 gives rise to (44)—(47). 


ias + iar = fate + jam (44) 
fog + igr = “gga + İqm (45) 
İidfe = igs + tar — jam (46) 
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ids Rs Lis Oe Lisigs Lis Oe Ltrigr R; A 


Figure 5. D-axis equivalent circuit of a DFIG with single core resistance 
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igs 


Figure 6. Q-axis equivalent circuit of a DFIG with single core resistance 


Sets of current differential equations applied in the DFIG simulations are presented in (48) to (53). 


as = i (Vas — Rgigs + WeLrsigs — Reeiare) (48) 
os = = (Vas = Reigs — WeLrsias — Reeigte) (49) 
“ee = = (Var — Relay — WeLmiqm + Webprigr — Reiare ) (50) 
= = r (Var — Rrigr — WeLmidm + WeLrriar — Reeigte) GD 
Sam = = (Rrelare + Welmigm) (52) 
Tan = (Relate — Welmiam) (53) 


Similarly, derivatives of (46) and (47) for a change in the core current gave rise to (54) and (55). 
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dt dt + dt dt S ) 
digte _ digs " digr _ digm (55) 
dt dt dt dt 


Power loss minimization assessment of a doubly fed induction generator ... (Crescent Onyebuchi Omeje) 


312 m) ISSN: 2252-8792 


The estimated stator active and reactive powers are expressed in (56). 


P, = = [Vas igs + Vas r 


oe (56) 
Qs = 3 [Vasias E Vasigs| 
Also, the estimated rotor active and reactive powers are expressed in (57). 
3 : : 
P: = = [Var lar + Vor | (57) 
3 ; A 
Q = z [Variar = Varigr| 
The mechanical dynamics of the DFIG based on speed and torque are presented in (58)—(60). 
3 Se i dk 
Tem = 2 x Lm(iamiqr M iqmiar) (58) 
dor 1 
dan aa (59) 
1 
Or = Jz x (Te Te Bwr) (60) 


The efficiency of the machine in terms of power output and total losses of the DFIG is expressed in (61). 


71 E Power Output = TemX@r i x 100% (61) 


ra Tot 
Power Output+Total Losses (TemX@r) + Pp ote 


5. RESULTS AND DISCUSSION 

The simulation parameters are presented in Tables 1 and 2. A plot of the derived flux linkage against 
the d-axis rotor current at varying core resistance is shown in Figure 7(a). It is observed that the flux linkage 
increases in linear proportion with the d-axis current with a corresponding rise in the slope and attains a 
Rfe 
10. 
total core loss was observed. The deviation in the exponential curve with higher value of core loss (900 kW) 


higher value with a reduced equivalent core-resistance (Req = ). In Figure 7(b) an exponential rise in the 


: R RO a i . f : 
is more pronounced when Reg = a which implies that more field current and circulating flux is produced at 


a lower core-resistance. In Figure 8(a), a plot of the total resistive loss against current was presented at a 
varying rotor resistance. It is also observed that as the resistance value is increased to R} = 40 x 2.250, the 
total resistive loss obtained peaked at 2000kW at a current value of 3.25A as against 3.5A got when the rotor 
resistance was reduced which accounted for the higher current flow through the path of least resistance. In 
Figure 8(b), it is observed that the DFIG total power loss rose to a value of 3250 kW at a reduced core 


: R : 4 : bos ; 
resistance value of Req = an and higher rotor resistance value of R, = 40 x 2.250 with a rise in slope. This 


indicates that at reduced core resistance, more field current is drawn since more paralleled resistance is 
needed to achieve this process whereas with an increased rotor resistance, more power is dissipated across 
the rotor terminal. The dq-axes voltages for the grid side converter (GSC) and rotor side converter (RSC) are 
presented in Figures 9(a) and 9(b). It is observed that the voltage magnitude is maintained at 450V during the 
simulation period thereby ensuring system stability in supply. In Figures 10(a) and 10(b), dq-axes rotor currents 
are shown at variable core resistance values of Rfe = 0.750 and 0.250 Similarly, in Figures 11(a) and 11(b), 
the plots obtained for the dq-axes core current showed that transient disturbance has more impacts on the 
DFIG core winding at a reduced core resistance value of Rfe = 0.250 and this accounted for its low 
efficiency value. In Figures 12(a) and 12(b), a sharp deviation in the settling time of the dq-axes magnetizing 
currents were observed at different core resistance values. The plots for the motor speed and electromechanical 
torques for different core resistance values are shown in Figures 13(a) and 13(b). It is observed that the transient 
response in speed and torque is reduced during the period of transient disturbance for Rfe = 0.750. but more 
oscillatory transient response was obtained after the disturbance was cleared. The power losses and power 
outputs obtained during a transient disturbance at core resistance values of Rfe = 0.750. and 0.250. are shown 
in Figures 14(a) and 14(b). It is shown that more losses were obtained in Figure 14(b) with a value of 
757.325 kW than in Figure 14(a) that dissipated 202.274 kW. This accounted for the large variation in their 
efficiency value of 83.45% and 41.21%. 
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In Figures 15(a) and 15(b), the stator active and reactive powers are presented. It is observed that a 


rapid oscillation was obtained during a transient disturbance which is more pronounced when the core 
resistance is kept at Rfe = 0.750. Conversely, in Figures 16(a) and 16(b), during a transient disturbance a 
rapid transient oscillation was observed at Rfe = 0.250 for rotor active and reactive power which implies 
that more current is drawn when Rfe = 0.250. 


A Plot of Flux Linkage against Current at varying Core Resistance 


1.8 


Table 1. Simulation parameters of the turbine blade model 


Parameters 


Values 


Blade radius (m) 


Maximum power coefficient (Cp) 


Optimal tip speed ratio 

Cut in speed (m/s) 

Rated wind speed (m/s) 
Minimum Pitch angle (degree) 
Maximum Pitch angle (degree) 


Table 2. Simulation parameters used for the DFIG dynamic model 


Parameters Values 
Rated power (MW) 15 
Rated mechanical power (MW) 10 
Supply Voltage (V) 400 
Supply Frequency (Hz) 50 
Stator resistance (pu) 0.25 
Rotor resistance (pu) 0.45 
Stator leakage inductance (pu) 0.0877 
Rotor leakage inductance (pu) 0.077 
Mutual inductance (pu) 0.955 
Core resistance of modified circuit (Q) 2.25 
Core resistance of simplified circuit (pu) 0.25 and 0.75 
Filter resistance (pu) 0.03 
Filter inductance (pu) 0.25 
Converter resistance (pu) 0.001 
Motor inertia Kg-M°) 0.045 
Coefficient of Viscosity (NMS) 0.00006 
Pole pair 2 
Load torque (NM) 20 


1.65 


0.8 F 


Minimized Flux Linkage (V.S) 


1000 
— Req = Rfe/10 900 f 
Req = Rfe/8 800 L 

—— Req = Rfe/4 
—— Req = Rfe/2 700 + 
600 + 
500 + 


A Plot of DFIG Core Loss against minimized Flux linkage (V.S) 
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Figure 7. A plot of (a) flux linkage (V.S) against current (A) and 
(b) total core loss (W) against flux linkage (V.S) 
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Figure 8. A plot of (a) total resistive loss (W) against current (A) and 
(b) DFIG total loss (W) against flux linkage (V.S) 


A Plot of d-axis Stator Voltage (V) against Time (Sec.) A Plot of d-axis Rotor Voltage (V) against Time (Sec.) 
r T T T T T 


500 500 
=o = 
= 5 
7 > 
-500 i 4 L 1 i 1 1 
0 006 O01 01 02 0235 03 035 04 -500 d i l j i £ 
Time( Sec.) i) 0.05 0.1 0.15 02 0.25 03 0.35 04 
A Plot of q-axis Stator Voltage (V) against Time (Sec.) E Time (Sec.) A 
500 T T T T T T T s00 A Plot of q-axis Rotor Voltage (V) against Time (Sec.) 
r 1 T T - 7 7 
= 0 = 
o = 0 
> ica 
> 
-500 1 A A A f a 
y 1 J H x p E 
° ae n uP 1 wa 3 s 4 0 0.05 0.1 0.15 02 0.25 03 0.35 04 
Time(Sec.) Time (Sec.) 
(a) (b) 


Figure 9. A plot of (a) dq-axes GSC voltage (V) against time (S) and 
(b) dq-axes RSC voltage (V) against time (S) 
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Figure 10. A plot of (a) dq-axes rotor current (A) for Rfe = 0.75Q and 
(b) dq-axes rotor current for Rfe = 0.250 
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A Plot of d-axis Core Current Idfe (A) against Time (Sec.’ A Plot of d-axis Core Current Idfe (A) against Time (Sec.) 
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Figure 11. The plots obtained for (a) dq-axes core current (A) for Rfe = 0.75Q. and 
(b) dq-axes core currents for Rfe = 0.250 
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Figure 12. A sharp deviation in the settling time of (a) dq-axes magnetizing current (A) for Rfe = 0.75Q and 
(b) dq-axes magnetizing currents for Rfe = 0.250 
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Figure 13. A plot of (a) speed (Rad/Sec) and torque (Nm) for Rfe = 0.75Q and (b) for Rfe = 0.250 
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A Plot of DFIG Power Output (kW) against Time (Sec.) 
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A Plot of DFIG Power Loss (kW) against Time (Sec.) 
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Figure 14. A plot of power loss and power output (kW) for (a) Rfe=0.75Q and (b) Rfe = 0.250 
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A Plot of Stator Active Power (kW) against Time (Sec.) 
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A Plot of Stator Reactive Power (kVAR) against Time (Sec.) 
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Figure 15. A plot of stator active and reactive power for (a) Rfe = 0.75Q and (b) Rfe = 0.25Q 


A Plot of DFIG Rotor Active Power (kW) against Time (Sec.) 
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A Plot of DFIG Rotor Active Power (kW) against Time (Sec.) 
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A Plot of DFIG Rotor Reactive Power (kVAR) against Time (Sec. 
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Figure 16. A plot of rotor active and reactive power for (a) Rfe = 0.75Q. and (b) Rfe = 0.25Q 
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6. CONCLUSION 

In this paper, power loss minimization of a DFIG has been presented using a modified DFIG 
equivalent circuit with multi-core resistances connected in parallel. The impact of core arrangements on 
dynamic power loss inside a DFIG was analyzed. The loss minimization model of the DFIG was developed 
incorporating iron, copper loss components, RL-filter losses, Frictional losses and power electronics 
converter losses with a minimum flux linkage to minimize the overall total power losses. Simulation results 
showed that losses are minimized when the equivalent core resistances are connected in parallel with 
minimum permissible current flow. This aids in regulating the magnitude of the reactive power between the 
wind turbine and the DFIG. An unconstrained optimization carried out showed that the DFIG power loss 
model analyzed using the Hessian matrix is positive definite and therefore is controllable. The results 
obtained during a transient disturbance indicated that at different core resistance values of Rfe = 0.750 and 
0.250. different efficiency values of 83.45% and 41.20% were obtained due to the variation in the core 
currents and core losses. It can be concluded from the simulation results that with the reduced core resistance 
and increased flux linkage more current flows through the core and more losses are dissipated with a 
consequent local saturation of the core which may be unavoidable. 
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